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Both hydrogen sulfide (H 2 S) and methyl mercaptan (CH 3 SH) are frequently detected in large amounts in 
malodorous mouth air. We investigated the bacterial composition of saliva of 30 subjects with severe oral 
malodor exhibiting extreme CH 3 SH/H 2 S ratios (high H 2 S but low CH 3 SH concentrations, n = 14; high 
CH 3 SH but low H 2 S concentrations, n = 16) and 13 subjects without malodor, using barcoded 
pyrosequencing analysis of the 16S rRNA gene. Phylogenetic community analysis with the UniFrac distance 
metric revealed a distinct bacterial community structure in each malodor group. The H 2 S group showed 
higher proportions of the genera Neisseria, Fusobacterium, Porphyromonas and SRI than the other two 
groups, whereas the CH 3 SH group had higher proportions of the genera Prevotella, Veillonella, Atopobium, 
Megasphaera, and Selenomonas. Our results suggested that distinct bacterial populations in the oral 
microbiota are involved in production of high levels of H 2 S and CH 3 SH in the oral cavity. 



Oral malodor is a common but distressing condition that may affect interpersonal social communication 1 . 
Most of the pathologies are found within the mouth and originate mainly from microbial metabolism 1 . 
The intraoral surfaces are colonized by a diverse array of bacterial species, many of which, especially 
anaerobes, have an ability to degrade substrates into malodorous compounds 2,3 . Microbial overgrowth due to 
poor oral hygiene results in the development of oral malodor. 

Molecular approaches using the 16S rRNA gene have allowed comprehensive surveys of complex bacterial 
communities, and the tongue microbiota structure of subjects with severe oral malodor has been determined 4 " 6 . 
These studies indicated that the microbiota of subjects with oral malodor shows a greater microbial diversity than 
that of no-odor subjects, and various microorganisms have been implicated in oral malodor. Our previous study 
of 240 subjects using terminal restriction fragment length polymorphism (T-RFLP) analysis confirmed that 
Gram-positive saccharolytic species, including Streptococcus, are more predominant in the microbiota of subjects 
with no odor than in that of malodor patients. Furthermore, our data demonstrated statistically that the global 
composition of the oral microbiota is correlated with the severity of malodor 7 . However, no single specific 
bacterial infection has been definitively associated with oral malodor. Key bacterial species to be targeted by 
therapy remain poorly defined. Some researchers have proposed the "non-specific theory" for the microbial 
etiology of oral malodor, which suggests that there are multiple causal agents and that many groups can substitute 
for others 8 . 

Unpleasant oral odor results from a mixture of various molecules, such as volatile sulfur compounds (VSC), 
short-chain fatty acids, polyamines, and indoles 8 . The oral malodor- related microbiota is presumed to be a 
complex consortium of bacterial members involved in production of these compounds. Discrimination of the 
various production systems might help to identify those microbes important in generation of oral malodor. 

Among the various compounds causing oral malodor, two VSCs, hydrogen sulfide (H 2 S) and methyl mer- 
captan (CH 3 SH), are known to be associated with oral malodor, and their concentrations in mouth air are often 
used for clinical diagnosis. In this study, we aimed to distinguish and identify the oral microbiota structures 
associated with high H 2 S and CH 3 SH production. Although malodorous mouth air contains both H 2 S and 
CH 3 SH in most cases, some patients exhibit extreme CH 3 SH/H 2 S ratios. We enrolled two groups of severe 
halitosis patients with distinct VSC profiles of high H 2 S but low CH 3 SH concentrations and high CH 3 SH but 
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Figure 1 | Hydrogen sulfide (H 2 S) and methyl mercaptan (CH 3 SH) 
distribution in mouth air of 240 subjects in our previous study 7 . Subjects 
with low VSC in mouth air (•, H 2 S < 0.075 ppm and CH 3 SH = 0 ppm, 
n = 13), those with high H 2 S but low CH 3 SH concentrations (V, H 2 S > 
1 ppm and CH 3 SH/H 2 S < 0.6, n = 14), and those with high CH 3 SH but 
low H 2 S concentrations (A, CH 3 SH > 0.5 ppm and CH 3 SH/H 2 S > 1.0, 
n = 16) were enrolled in this study. The other subjects shown as gray circles 
were excluded from the present study. 

low H 2 S concentrations, as well as a group of control subjects with no 
oral malodor. The bacterial composition of saliva was evaluated 
using barcoded pyrosequencing analysis of the 16S rRNA gene. 

Results 

We investigated the salivary bacterial populations of 43 subjects 
whose mouth air showed three different VSC profiles: high H 2 S 
but low CH 3 SH concentration (H 2 S group, n = 14), high CH 3 SH 
but low H 2 S concentration (CH 3 SH group, n = 16), and no VSC (no- 
odor group, n = 13) (Figure 1). The general and clinical parameters 



Table 1 General 


and oral condition of 43 subjects 

No odor group H2S group 
(n = 1 3) (n = 14) 


CH3SH group 
(n= 16) 




28±11 


47±15 b 


57±1 3 b 


oex (remale [/o]) 


6(46.1) 


5 (35.7) 


O KA 0\ 
V po.zj 


IN umber or teeth 


27±4 


26±4 


0 R ■+- Q 


INUIIlUci vJI 


1.3±4.6 


1.4±3.1 


W.O — \J.vJ 


decayed teeth a 








Tongue coating 


1.2±0.5 


2.0±0.8 b 


2.1 ±0.8 b 


score a 








Mean pocket 


2.8±0.2 


3.4±0.9 


3.4±0.5 b 


depth a (mm) 








Amount of saliva 0 


5.7±2.5 


7.3±3.9 


6.3±2.5 


(ml/5 min) 








"Significant difference between groups were evaluated by Steel-Dwass multiple comparison test. 
Significantly higher than no odor group (P < 0.05). 
Significant difference were evaluated by Fisher's exact test. 



of the study populations are given in Table 1 . Age and tongue coating 
scores of the two malodor groups were significantly higher than those 
of the no-odor group, whereas no significant difference was observed 
between the two malodor groups. In addition, the mean periodontal 
pocket depth of the CH 3 SH group was significantly greater than that 
of the no-odor group. 

Whole-genomic DNA was extracted from each saliva sample and 
the variable regions (VI -V2) of the bacterial 16S rRNA gene were 
PCR amplified using barcoded universal primers. Pyrosequencing 
using the Roche 454 FLX instrument produced a dataset consisting 
of 231,266 high-quality reads with an average length of 344 ± 20 bp 
(Table 2). Of these, 196,943 were matched with 575 oral bacterial 
sequences deposited in the Forsyth Institute Human Oral Micro - 
biome (HOMD) database 9 , and these sequences were subsequently 
assigned to 429 oral taxon (OT) defined by HOMD 10 . The remaining 
34,323 sequences were assigned to 2,639 operational taxonomic units 
(OTU) using a 97% pairwise-identity cutoff. 

Phylogenetic community analysis with the UniFrac distance metric 
revealed a distinct overall bacterial community composition in each 
of the three groups (Figure 2). These differences were confirmed 
statistically in both unweighted (qualitative) and weighted (quantitat- 
ive) versions of this metric using permutational multivariate analysis 



Table 2 | 


Summary of pyrosequencing analysis 












No odor group 




H2S group 




CH3SH group 




(n= 13) 




(n= 14) 




(n= 16) 


Subject ID 


No. of Reads 


Subject ID 


No. of Reads 


Subject ID 


No. of Reads 


Nl 


4,249 


HI 


8,661 


Ml 


5,825 


N2 


6,858 


H2 


7,667 


M2 


3,720 


N3 


7,044 


H3 


5,278 


M3 


2,714 


N4 


5,552 


H4 


6,310 


M4 


6,748 


N5 


6,273 


H5 


3,077 


M5 


6,412 


N6 


3,912 


H6 


1,655 


M6 


2,799 


N7 


7,076 


H7 


4,81 1 


M7 


6,033 


N8 


6,605 


H8 


5,943 


M8 


5,564 


N9 


4,673 


H9 


7,218 


M9 


4,987 


N10 


3,591 


H10 


7,087 


M10 


6,677 


Nl 1 


5,709 


HI 1 


4,336 


Ml 1 


3,649 


N12 


4,417 


H12 


3,745 


M12 


3,904 


N13 


5,553 


H13 


7,031 


M13 


5,689 






H14 


5,590 


M14 


5,224 










M15 


3,792 










M16 


7,608 


Average 


5,500± 1,230 




5,600± 1,938 




5,084± 1,490 
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Figure 2 | Principal coordinate analysis (PCoA) plot showing similarity of composition of 43 salivary microbiotas. Plots were generated using 
unweighted (A) and weighted (B) versions of the Unifrac distance metric. The two components explained 15.9 and 71.7% of the variance, respectively. 

of variance (perMANOVA) and analysis of similarities (ANOSIM) metric (R = 0.25), suggesting that the microbiotas of the three groups 
analyses (P < 0.001, each). The discrimination was stronger with the of subjects had greater differences in community structure than in 
weighted UniFrac (ANOSIM statistic R = 0.47) than the unweighted community membership. 
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Figure 3 | The relative abundances of phyla in the salivary microbiota of 43 subjects (mean ± SD). Significant differences between groups were 
evaluated by Steel-Dwass test. **P < 0.01, *P < 0.05. 
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The vast majority of the sequences were assigned to five bacterial 
phyla (Actinobacteria, Bacteroidetes, Fusobacteria, Firmicutes, and 
Proteobacteria). TM7, SRI, Spirochaetes, Tenericutes, Cyanobacteria, 
and Synergistetes were also identified but in much lower proportions. 
Significant differences were observed in the relative proportions of 
the five major phyla and SRI among the three groups (Figure 3). 

The sequences represented 82 different genera (constituting at 
least 97% of the microbiota of each subject) and 35 upper-level taxa. 
Of the 82 genera, 58, 66, and 68 were detected in the no-odor, H 2 S, 
and CH 3 SH groups, respectively. The overall composition of the 
microbiota was qualitatively fairly uniform, and 49 genera were 
detected in each group. Eleven genera (Streptococcus, Prevotella, 
Rothia, Actinomyces, Granulicatella, Neisseria, Terrahaemophilus, 
Veillonella, Gemella, Fusobacterium, and Haemophilus) were found 
in all subjects (Figure 4) and constituted 80 to 99% of the microbiota 
in each individual. On the other hand, the relative abundance of each 
genus varied among the three groups. Significant differences between 
the three groups were observed in 19 bacterial genera (Figure 5). 
Compared to the no- odor group, the microbiota of both malo- 
dor groups was colonized by significantly lower proportions of 
Streptococcus and Granulicatella and higher proportions of Leptotri- 
chia, Peptostreptococcus, Eubacterium, and Fusobacterium. On the 
other hand, the genera Fusobacterium, Neisseria, Porphyromonas, 
and SRI genera incertae sedis were significantly more predominant 
in the H 2 S group than in the other two groups. Conversely, the genera 
Prevotella, Veillonella, Atopobium, Megasphaera, and Selenomonas 
were significantly more predominant in the CH 3 SH group than 
in the other two groups. Significant differences were also observed 
in Rothia (no-odor group > H 2 S group), Parvimonas (H 2 S group > 



No odor group 




H 2 S group 




Figure 4 | The mean abundances of bacterial genera in the no-odor, H 2 S, 
are shown. 



no-odor group), Peptococcus (H 2 S group > no-odor group), Solo- 
bacterium (CH 3 SH group > no-odor group), and Dialister (CH 3 SH 
group > no-odor group). 

Of 429 OT and 2,639 OTU, 77, including well-known VSC pro- 
ducers such as Porphyromonas gingivalis (OT 619), P. endodontalis 
(OT 273), and Fusobacterium nucleatum (OT 200 and OT 202), were 
found in significantly higher proportions in both malodor groups 
than in the no-odor group (Fig. 6 and Supplementary Table SI). In 
contrast, 15 of these were significantly more predominant in the 
microbiota of the H 2 S group compared to the CH 3 SH group; con- 
versely, 17 were significantly more predominant in the microbiota of 
the CH 3 SH group compared to the H 2 S group (Figure 6). 

Discussion 

This study revealed the bacterial community structure in the saliva of 
two groups of severe malodor patients with different CH 3 SH/H 2 S 
ratios. The microbiota of the subjects in both malodor groups was 
characterized by lower proportions of Streptococcus and Granuli- 
catella and higher proportions of Leptotrichia, Peptostreptococcus, 
Eubacterium, and Fusobacterium compared with the no-odor group. 
On the other hand, the genera Neisseria, Fusobacterium, Porphyro- 
monas, and SRI and the genera Prevotella, Veillonella, Atopobium, 
Megasphaera, and Selenomonas were detected in higher proportions 
in the H 2 S and CH 3 SH groups, respectively (Figure 5). Many oral 
bacteria can produce both H 2 S and CH 3 SH 2 , and their concen- 
trations in mouth air are generally highly correlated each other. 
Therefore, previous studies explored the oral malodor-associated 
microbiota without focusing on the VSC distribution. However, 
our results clearly indicate that the bacterial community structures 
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Figure 5 | The relative abundances of 19 bacterial genera that differed significantly between groups. Statistical differences were evaluated using a Steel- 
Dwass multiple comparison test (P < 0.05). Significant difference between the CH 3 SH and no-odor groups. "Significant difference between the H 2 S and 
no-odor groups. Significant difference between the CH 3 SH and H 2 S groups. 



of the two malodor groups were phylogenetically different 
(Figure 2), suggesting that distinct bacterial populations are 
involved in production of high levels of each compound in the oral 
cavity. 

Whereas potential H 2 S producers such as Fusobacterium and 
Porphyromonas were significantly more predominant in the H 2 S 
group, all characteristic bacterial species, especially in the CH 3 SH 
group were not necessarily heavy producers of each compound. 
According to Persson et al. 2 , many Prevotella species, including P. 
melaninogenica, have no, or only a low, capacity to produce CH 3 SH 
from serum or methionine. Veillonella and Selenomonas species 
degrade neither serum nor methionine into significant amounts of 
CH 3 SH, although they are potent H 2 S producers. To our knowledge, 
no report has demonstrated that Atopobium and Megasphaera spe- 
cies have a high capability to generate CH 3 SH. VSC production in the 
oral cavity is considered to result from proteolysis of serum, food 
debris and other substrates, followed by subsequent breakdown of 
sulfur- containing amino acids 11 . VSC production is likely to be a 
complex process that involves many bacterial species including both 
VSC producing and non-producing strains. Sterer and Rosenberg 
indicated that even saccharolytic species such as Streptococcus sali- 
varius promote VSC production through supply of protein resulting 



from glycoprotein putrefaction 11 . Well-known producers of both 
VSCs, such as Porphyromonas gingivalis, P. endodontalis, and Fuso- 
bacteirum nucleatum, were detected in higher proportions in both 
malodor groups compared with the no-odor group (Figure 6), 
although they represented only small minorities of the total micro- 
biota (Supplementary Table SI). Their H 2 S or CH 3 SH productivity 
might be accelerated by synergistic interactions with the surrounding 
bacterial populations. There is also the possibility that these bacteria 
are non-participating observers that prefer microbial environments 
in which production of large amounts of H 2 S or CH 3 SH occurs. Co- 
cultivation in vitro would assist in identification of the key bacterial 
combinations. In addition, the dominance of SRI observed in the 
H 2 S group (Figures 5 & 6) is noteworthy, considering that the eco- 
logical role of this as -yet- uncultured phylum is not well under- 
stood 12 . Development of a method of cultivating this bacterium 
would clarify its role in high H 2 S production. 

Our previous study using T-RFLP classified the salivary microbiota 
of 240 subjects into four patterns. Two of these (clusters II and IV) were 
equally implicated in being associated with oral malodor 7 . The char- 
acteristics of clusters II and IV, which were predicted from their T- 
RFLP profiles, corresponded well to the bacterial community structures 
of the CH 3 SH and H 2 S groups, respectively, in the present study. 
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Figure 6 | Relative abundance distribution of 77 oral taxon (OT) or operational taxonomic units (OTU) that were significantly more abundant in the 
malodor groups (n = 30) than in the no-odor group (n = 13), evaluated by Wilcoxon rank-sum test (P < 0.05). To show the distribution of the OT or 
OTU with lower abundance, the relative abundances of each OT or OTU were normalized to have a mean of 0 and standard deviation of 1 (z-score 
normalization) and represented as the color intensity of each grid (blue, low abundance; red, high abundance). The OT and OTU were ordered according 
to the Z statistic (shown as a bar plot on the left) in a Wilcoxon rank-sum test comparing the CH 3 SH (n = 16) and H 2 S (n = 14) groups. The direction of 
horizontal bars indicates greater predominance in the H 2 S than the CH 3 SH group (left) and in the CH 3 SH than the H 2 S group (right). The red and blue 
bars show significant differences between the two malodor groups. 



However, their association with CH 3 SH and H 2 S was not examined in 
the previous study. Therefore, we reanalyzed the data regarding CH 3 SH 
and H 2 S in each cluster. The results indicated that the adjusted odds 



ratio for increased CH 3 SH levels was higher in cluster II than cluster IV, 
and conversely, that the odds ratio for increased levels of H 2 S were 
higher in cluster IV than cluster II (data not shown). These findings are 
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consistent with the relationship between high H 2 S and CH 3 SH produc- 
tion and the microbiota structures characterized in this study. 

Along with H 2 S and CH 3 SH, various other malodorous com- 
pounds are often present in malodorous mouth air, such as short- 
chain fatty acids, polyamines, and indoles 8 . The oral microbiota of 
the high VSC groups contained not only VSC-producing bacteria, 
but also those capable of producing other malodorous compounds. 
Therefore, further subtyping of high VSC-producing microbiota 
based on the yields of other malodorous compounds might facilitate 
discrimination and identification of bacterial members involved in 
their production, if appropriate clinical methods of evaluating these 
compounds were available. Elucidating the bacterial species involved 
in production of malodorous compounds other than VSC would 
help to further narrow the key bacterial species responsible for 
VSC production. 

Consistent with previous studies, we found no single bacterial 
species whose presence explains oral malodor. However, our results 
revealed that a large amount of H 2 S or CH 3 SH in the oral cavity is 
produced by characteristic oral microbiota structures, suggesting 
that suppression or elimination of specific bacterial populations 
may effectively reduce oral malodor. Although a methodology of 
altering the stable oral microbiota remains unknown, this may rep- 
resent a treatment option in addition to the conventional steadfast 
mouth- cleaning routine that focuses on reducing the total bacterial 
load in the oral cavity. 

Methods 

Study population. The subjects were a subgroup of the 240 subjects analyzed in our 
previous study (Takeshita et al. 2010), who visited the Oral Malodor Clinic at 
Fukuoka Dental College Medical and Dental Hospital in Japan. In this study, we 
selected 30 subjects with severe oral malodor showing two characteristic VSC profiles: 
high H 2 S but low CH 3 SH concentration (H 2 S > 1 ppm and CH 3 SH/H 2 S < 0.6, 
n = 14) and high CH 3 SH but low H 2 S concentration (CH 3 SH > 0.5 ppm and 
CH 3 SH/H 2 S > 1.0, n = 16) (Figure 1). We also enrolled 13 subjects without oral 
malodor (H 2 S < 0.075 ppm and CH 3 SH = 0 ppm) based on the H 2 S and CH 3 SH 
olfactory thresholds (0.15 and 0.05 ppm, respectively 13 ). Their general and clinical 
conditions were evaluated in the previous study 7 . The Ethics Committee of Fukuoka 
Dental College and Kyushu University Faculty of Dental Science approved the study 
design. 

Sample collection and pyrosequencing analysis. Collection of saliva samples and 
DNA extraction was performed in the previous study (Takeshita et al. 2010). The 16S 
rRNA genes of each sample were amplified using the following primers: 338R with the 
454 Life Sciences adaptor B sequence (5' -CCT ATC CCC TGT GTG CCT TGG CAG 
TCT CAG TGC TGC CTC CCG TAG GAG T-3') and 8F with the 454 Life Sciences 
adaptor A and subject- specific six-base barcode sequences (5' -CCA TCT CAT CCC 
TGC GTG TCT CCG ACT CAG NNN NNN AGA GTT TGA TYM TGG CTC 
AG-3'). PCR amplification was performed as described previously (Takeshita et al. 
2010). The amplicons were gel-purified using a Wizard SV Gel and PCR Clean-Up 
System (Promega, Madison, WI) according to the manufacturer's instructions. The 
DNA concentration and quality were assessed using a NanoDrop spectrophotometer 
(NanoDrop Technologies, Wilmington, DE), and equal amounts of DNA from the 43 
subjects were pooled. Pyrosequencing was conducted using a 454 Life Sciences 
Genome Sequencer FLX instrument (Roche, Basel, Switzerland) at Hokkaido System 
Science Co., Ltd. (Sapporo, Japan). 

Data analysis and taxonomy assignment. Sequences were excluded from the 
analysis using a script written in PHP if they were shorter than 240 bases, or had an 
average quality score <25, and subsequently removed using a script written in R if 
they did not include the correct primer sequence, had a homopolymer run >6 nt, or 
contained ambiguous characters. The remaining sequences were assigned to each 
subject by examining the six-base barcode sequence. Each sequence was compared 
using the BLAST algorithm with 1,647 sequences of the oral bacterial 16S rRNA gene 
(HOMD 16S rRNA RefSeq Extended Version 1.1) deposited in HOMD 9 , and 
assigned to the best BLAST hit with a 98.5% similarity value and minimum coverage 
of 97%. After excluding assigned sequences, similar sequences were clustered into 
operational taxonomic units using the complete-linkage clustering tool of the RDP 
pyrosequencing pipeline 14 at a distance cut off of 0.03, and the representative 
sequences of each cluster were selected using the Dereplicate request function. The 
oral bacterial sequences from HOMD, to which at least one sequence was assigned, 
and the representative sequences from each OTU were aligned using PyNAST 15 and 
the Greengenes database 16 using a minimum percent identity of 75%. Chimeras were 
removed from the representative set on the basis of identification as chimeric via 
Chimera Slayer 17 and verification that the putative chimera appeared in only one 
sample. After chimera elimination, a relaxed neighbor- joining tree was built using 



FastTree 18 . To determine the dissimilarity between any pair of bacterial communities, 
we used the UniFrac metric 19 calculated by Fast UniFrac 20 . The similarity relationship 
assessed using the unweighted and weighted UniFrac metric was represented in a 
principal coordinate analysis (PCoA) plot drawn by R. The taxonomy of 
representative sequences was determined using the RDP classifier with a minimum 
support threshold of 60% and the RDP taxonomic nomenclature (down to the genus 
level). 

Statistical analysis. A Steel-Dwass non-parametric multiple comparison test was 
performed to compare age, number of teeth, number of decayed teeth, tongue coating 
score, mean pocket depth, amount of saliva, and the relative abundances of bacterial 
phyla and genera. Fisher's exact test was conducted to look for differences by sex. 
Wilcoxon's rank- sum test was performed to compare the relative abundance of each 
OT or OTU. We used perMANOVA and ANOSIM to test for differences in bacterial 
community structure among groups of samples. Statistical significance was set at 
P < 0.05. All statistical analyses were conducted using the R software package, version 
2.11. 1 21 . 
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